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We report observations of superradiance for atoms trapped in the near field of a photonic crystal waveguide 
(PCW). By fabricating the PCW with a band edge near the Di transition of atomic cesium, strong interaction 
is achieved between trapped atoms and guided-mode photons. Following short-pulse excitation, we record the 
decay of guided-mode emission and find a superradiant emission rate scaling as fsR oc iV • Fid for average 
atom number 0.19 < N < 2.6 atoms, where Fid/Fo = 1.1 ± 0.1 is the peak single-atom radiative decay rate 
into the PCW guided mode and Fq is the Einstein-A coefficient for free space. These advances provide new 
tools for investigations of photon-mediated atom-atom interactions in the many-body regime. 


Interfacing light with atoms localized near nanophotonic 
structures has attracted increasing attention in recent years. 
Exemplary experimental platforms include nanofibers iiiia, 
photonic crystal cavities a and waveguides naa. Owing 
to their small optical loss and tight field confinement, these 
nanoscale dielectric devices are capable of mediating long- 
range atom-atom interactions using photons propagating in 
their guided modes. This new paradigm for strong interac¬ 
tion of atoms and optical photons offers new tools for scalable 
quantum networks 0, quantum phases of light and matter 
lllIHl, and quantum metrology |j9l. 

In particular, powerful capabilities for dispersion and 
modal engineering in nanoscopic photonic crystal waveguides 
(PCWs) provide opportunities beyond conventional settings 
in AMO physics within the new field of waveguide QED 
fflimiiiioHni- For example, the edge of a photonic band gap 
aligned near an atomic transition strongly enhances single¬ 
atom emission into the one-dimensional (ID) PCW due to a 
van-Hove singularity at the band edge (i.e., a ‘slow-light’ ef¬ 
fect GMSl). Because the Bloch function for a guided mode 
near the band edge approaches a standing-wave, symmet¬ 
ric optical excitations can be induced in an array of trapped 
atoms, resulting in superradiant emission csiini into the 
PCW. Superradiance has important applications for realizing 
quantum memories 1(1814^ . single photon sources ||^|24l, 
laser cooling by way of cooperative emission 125] Hi, and 
narrow linewidth lasers EU. Related cooperative effects 
are predicted in nano-photonic waveguides absent an exter¬ 
nal cavity ll28l . including atomic Bragg mirrors ll29l and self¬ 
organizing crystals of atoms and light ll30l43^ . 

Complimentary to superradiant emission is the collective 
Lamb shift induced by proximal atoms virtually exchanging 
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off-resonant photons Il3?l43^ . With the atomic transition fre¬ 
quency placed in a photonic band gap of a PCW, real photon 
emission is largely suppressed. Coherent atom-atom interac¬ 
tions then emerge as a dominant effect for QED with atoms in 
bandgap materials ETM). Both the strength and length scale 
of the interaction can be ‘engineered’ by suitable band shap¬ 
ing of the PCW, as well as dynamically controlled by external 
lasers HU] Ha. Exploration of many-body physics with tun¬ 
able and strong long-range atom-atom interactions are thereby 
enabled 114111421 . 

In this Letter, we present an important advance for the field 
of waveguide QED. We describe an experiment that cools, sta¬ 
bly traps, and interfaces multiple cold atoms along a quasi 
one-dimensional PCW. Through precise band edge alignment 
and guided-mode (GM) design, we achieve strong radiative 
coupling of one trapped atom and a GM of the PCW, such 
that the inferred single-atom emission rate into the GM is 
Fio/ro = 1.1 ± 0.1, where Fid is the peak single-atom ra¬ 
diative decay rate into the PCW guided mode and Fg is the 
Einstein-A coefficient for free space. With multiple atoms, we 
observe superradiant emission in both time and frequency do¬ 
mains with measurements of transient decay following pulsed 
excitation and steady-state transmission spectra, respectively. 
We infer cooperative, superradiant coupling with rate Fsr that 
scales with the mean atom number (V as FgR = r]N ■ Fid over 
the range 0.19 ^ N < 2.6 atoms, where 77 = 0.34 ± 0.06. 

Our experimental platform is based on trapped cesium 
atoms near a ID alligator photonic crystal waveguide 
(APCW) a El. The APCW is formed by two parallel SiN 
nanobeams separated by 238 nm with periodic corrugations 
at the outer edges (Fig. 1(a)). The APCW consists of 150 
identical unit cells with lattice constant a = 371nm (length 
L ~ 55.7/im) and is terminated at either end by 30 tapered 
cells for mode matching to parallel nanobeams without cor¬ 
rugation. Photons can be coupled into and out of the APCW 
from conventional cleaved-fibers at either end of the structure. 
Design principles, fabrication methods, and device character¬ 
ization of the APCW can be found in Refs. aiaiii. 

For the APCW used here, we align the band edge of 
the fundamental guided mode (electric field predominantly 
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transverse-electric (TE) polarized in the plane of the waveg¬ 
uide) near the cesium Di line at 894.6 nm, with a mode- 
matched TE input field Ein tuned around the 6 S 1 / 2 , F = 
3 —?► 6 Pi /2 5 F' = a transition. Near the band edge, the atom- 
photon coupling rate is significantly enhanced by the group in¬ 
dex Tig, as well as by reflections from the tapering regions that 
surround the APCW. Erom the measured transmission spec¬ 
trum of the device absent atoms, we estimate a group index 
rig ~ 11 and an intensity enhancement ~ 6 from the taper 
reflections El. 

To trap atoms along the APCW, we create tight optical po¬ 
tentials using the interference pattern of a side-illumination 
(SI) beam and its reflection from the surface of the APCW 
EEl. The polarization of the SI beam is aligned paral¬ 
lel to the x-axis of the ID waveguide to maximize the re¬ 
flected field. Eigure [TJb) shows the calculated near-fleld in¬ 
tensity distribution in the y-z plane ll45l . With a red-detuned 
SI beam, cold atoms can be localized to intensity maxima 
(e.g., positions Z-i,zi^Z 2 in Eig. [^b)). However, because 
of the exponential falloff of the GM intensity, only those 
atoms sufficiently close to the APCW can interact strongly 
with guided-mode photons of the input held Ein, Pig- 
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FIG. 1: Trapping and interfacing atoms with a ID photonic crys¬ 
tal waveguide, (a) A side-illumination (SI) beam is reflected from 
an ‘alligator’ photonic crystal waveguide (APCW) to form a dipole 
trap to localize atoms near the APCW (gray shaded structure). The 
red shaded region represents trapped atoms along the APCW. An in¬ 
cident field Ein excites the TE-like fundamental mode and thereby 
trapped atoms couple to this guided mode (GM). The transmitted 
tEin and reflected field rEi^ are recorded. The inset shows an SEM 
image of the APCW and corresponding single-atom coupling rate 
r ID along the x axis at the center of the gap (y = 0). (b) Normal¬ 
ized intensity cross section of the total intensity /tot resulting from 
the SI beam and its reflection, which form an optical dipole trap. 
Trap locations along the 2 axis at 1 / = 0 are marked by Zi. Masked 
gray areas represent the APCW. (c) The single-atom coupling rate 
into the TE guided mode riD(0, y, z) normalized to the free-space 
decay rate Fq for the cesium Di line. 



EIG. 2: Lifetime of trapped atoms near the APCW. (a) 1/e-lifetime 
of Tfa = 54±5 ms is determined using free-space absorption imaging 
of the trapped atom cloud, (b) 1 /e-lifetime of tgm = 28 ± 2 ms is 
observed from the normalized transmission T/Tq of resonant GM 
probe pulses. 

The trap site with the strongest atom-photon coupling is lo¬ 
cated at ( 2 / 1 , Zi) = ( 0 , 220 ) nm, closest to the center of the 
unit cell and Az ~ 120 nm from the plane of the upper sur¬ 
faces of the APCW. Other locations are calculated to have 
coupling to the fundamental TE-like mode less than 1% of 
that for site Zi (e.g., the sites at z-i, Z 2 have intensity ratios 
/( 2 _i)//( 2 i) = 0.01 ,/(z2)//(2i) = 0.005). 

Along the x axis of the APCW, the dipole trap U{x, 0, zi) is 
insensitive to the dielectric corrugation within a unit cell and 
is nearly uniform to within < 2 % around the central region 
of the APCW. By contrast, atom emission into the fundamen¬ 
tal TE-like mode is strongly modulated with riD(a;, 0, zi) ~ 
Eid cos^(fcx) due to the Bloch mode function near the band 
edge of the APCW {k « tt/o), as shown in the inset of 
Pig. [2 a). Thus, even for atoms uniformly distributed along 
the X axis of the trapping potential, only those close to the 
center of a unit cell can strongly couple to the guided mode, 
greatly facilitating phase-matched symmetric excitation of the 
atoms. In our experiment, we have chosen a 50 ym waist for 
the SI beam to provide weak confinement along the x axis, 
with atoms localized near the central region (Ax ~ ±10 pm) 
of the APCW for the estimated temperature ~ 50pK from a 
time-of-flight measurement in free space. The SI beam for 
dipole trapping is 220 GHz red-detuned with respect to the 
D 2 line and has a total power of 50 mW for all measurements 
reported. 

Cold atoms from a MOT that surrounds the APCW IS are 
loaded into the dipole trap during an optical molasses phase 
(~ 5 ms) and then optically pumped to 6 S'i/ 2 , F = 3 for ^1 
ms. Atoms are held in the dipole trap for time fhoid relative to 
the end of the loading sequence, and then free-space absorp¬ 
tion imaging is initiated over the interval (fhoid, fhoid ± Atm) 
with Atm — 0.2 ms. We introduce the measured time 
tm = ihoid ± Afni/2, Centered in the measurement window. 
As shown in Pig. |^a), we measure a trap lifetime Tfg = 54±5 
ms and And a peak density po ~ 2 x 10^^ cm“^ near the 
APCW. The atom density p near the APCW can be adjusted 
over a wide range 0.06 ^ p/po ^ 1 by varying the duration 
of the MOT loading cycle while keeping all other procedures 
identical. 

To determine the lifetime for trapped atoms near the 
APCW, we again hold atoms for fhoid, and then launch Em 
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as a resonant GM probe in measurement interval ± Atm/2 
with Atm = 5 ms. From the recorded transmitted signals, we 
compute T/Tq, where Tq is the transmission without atoms. 
During the probe period, we also apply free-space repump 
beams, tuned to the D 2 , QS 1 / 2 , F — 4 ^ QP 3 / 2 J F' = 4 
resonance, to remove population in the 6 S^, F = 4, since 
the probe excites an open transition. Fig. I^b) shows T/Tq 
gradually recovering to T/Tq = 1 as tm increases, with a ht 
to the data giving a 1 /e—time of tgm = 28 ± 2 ms ll44l . tgm 
is consistently shorter than Tfg from free-space imaging, which 
might be attributed to increased heating from the stronger light 
intensity near the APCW, the effect of surface potentials, or 
outgassing from the silicon chip and structures that support 
the APCW. These contributions are being investigated in more 
detail. 

Our principal investigation of superradiance involves ob¬ 
servation of the transient decay of emission from an array of 
atoms trapped along the APCW. For a collection of W > 1 
atoms, superradiance is heralded by a total decay rate Ftot = 
Tsk + r^ot that is enhanced beyond the total rate of decay 
for one atom F^^j = Fid -f F'. Fsr is the W-dependent su¬ 
perradiant rate operationally determined from Ftot and Fj^j. 
Here, F' is the radiative decay rate into all channels other 
than the TE-like GM of the APCW. We numerically evalu¬ 
ate F'/Fq « 1.1 for an atom at the trap site zi in Fig. 
along the APCW, with Fq the free-space decay rate for the Di 
transition lfT4ll46l . Cooperative level shifts |iFdd| ^ Fid are 
neglected for the current configuration of our experiment ll44l . 

We record the temporal profiles of atomic emission into 
the fundamental TE-like GM following short-pulse (~10 ns 
EWHM), resonant excitations via Em- To ensure small popu¬ 
lation in the excited state, we choose a pulse intensity well be¬ 
low the saturation intensity (///sat < 0.1). After a time fhoid 
the excitation cycle is repeated every 500 ns for Atm = 6 ms, 
and detection events are accumulated for the reflected inten¬ 
sity jrE'inp by an avalanche photodiode (APD). We consider 
decay curves of GM emission at 15 ns < fe < 70 ns after the 
center of the excitation pulse (i.e., after the excitation pulse 
is sufficiently extinguished, te > 15 ns, and while the back¬ 
ground counts are negligible compared to the atomic emis¬ 
sion, fe ^ 70 ns HU). The total decay rate ftot is extracted 
by simple exponential fits as shown in the inset of Pig. |^a). 
The deviation from the exponential fit at te > 60 ns is due 
to the spatially varying coupling rate Fid cos^{kx), which is 
captured by a detailed model discussed later HD. 

Enhanced total decay rate with increasing atom number is 
clearly evidenced in Pig. [^a), where the atom number can 
be adjusted by varying trap hold time fhoid prior to the mea¬ 
surement. At the shortest measurement time = 3 ms with 
^hoid = 0 ms (i.e., the maximum number of trapped atoms), 
the measured total decay rate is largest at ftot/ro ~ 2.9. 
At tm = 63 ms much longer than the trap lifetime tgm = 
28 ± 2 ms, the total decay rate settles to ftot/ro ~ 2.0. This 
asymptotic behavior suggests that f tot at long hold time cor¬ 
responds to the single-atom decay rate ftot- 




FIG. 3: Decay rate and atom number dependence (a) Fitted total de¬ 
cay rate ftot normalized with free-space decay rate Fq (circles) as a 
function of measurement time tm- The solid line is a simple exponen¬ 
tial fit to determine the superradiant decay rate f sr/Fq = 1.1 ± 0.1 
and the single-atom decay rate ft/t/fo = 2.0 ± 0.1 with tsr = 
17 ± 3 ms. The inset shows the temporal profiles of normalized 
guided-mode emission /p//po (circles) with Ipo the peak emission. 
Exponential fits (solid curves): = 3 ms (red), 13 ms (green), and 

63 ms (blue). The black dashed curve shows a exponential decay 
with free-space decay rate Fq. (b) Fitted total decay rate ftot nor¬ 
malized with Fo as a function of mean number of trapped atoms N 
from a detailed model (44). We adjust N by changing the trap hold 
time (red circles) or atom loading time (blue circles). The black line 
is a linear fit to the combined data sets, giving Fsr = p • A • Fid 
with 77 = 0.34 ± 0.06. 

To determine quantitatively the superradiant and single¬ 
atom emission rates from our measurements of decaying GM 
emission, we present two different analyses that yield consis¬ 
tent results. First is a simple and intuitive analysis applied 
to Fig. 3(a) in which we employ an empirical exponential 
fit, Ftot(fm) = fsRe“*“AsR _|_ f^^^, with the superradiant 
fsR, single-atom f^ot, tsr characterizing decay of su¬ 
perradiance due to the atom loss. The fit yields the maximum 
superradiant rate Fsr/Fo = 1.1 ± 0.1 with tsr = 17 ± 3 
ms, and a reasonable correspondence to the measured de¬ 
cay rates ftot, as shown by the red curve in Fig. [^a). The 
asymptote f^Qj/Fo = 2.0 ± 0.1 gives the total single-atom 
decay rate. With F'/Fq ~ 1.1 determined numerically for an 
atom at trap site zi along the APCW (Fig. [^b)), we deduce 
Fid/Fq = 0.9 ± 0.1 for the single-atom decay rate into the 
GM of the APCW. 

To substantiate this simple emphirical model, our second 
analysis is a detailed number treatment based upon transfer 
matrix calculations m. Decay curves are generated for a 















4 


fixed number of atoms N distributed randomly along along 
the x-axis of the APCW with uniform probability density but 
with spatially varying spatially varying coupling Fid (a;) — 
Fid cos^{kx). These A^-dependent, spatially-averaged decay 
curves are further averaged over a Poisson distribution with 
mean atom number N, capturing the variation of atom num¬ 
ber N as we repeat experiments for data accumulation. Fitting 
to this model, we extract Fid/Fq = 1.1 ± 0.1 for measure¬ 
ments at long hold time (e.g., at = 63ms in Fig. [^a)). 
Since the intensity of the fluorescence from a single atom is 
spatially modulated by cos^(fca;), only an atom near the cen¬ 
ter of unit cell can strongly couple to the GM, resulting in the 
small difference between averaged Fid and peak Fid. Also, 
the decay curve for GM emission at fm = 3 ms can be well fit¬ 
ted with TV = 2.6±0.3atoms ll44l . The red points inFig.|^(b) 
display the total decay rate ftot as a function of N extracted 
from fits of the transfer matrix model to the measured decay 
curves, which clearly shows that supetTadiance emission rate 
is proportional to N. 

The value Fid/Fq = 1.1 ± 0.1 from our measurements 
agrees reasonably well with the theoretical value Fid/Fq « 
1.2 determined by FDTD calculations ll43l l44l . despite sev¬ 
eral uncertainties (e.g., locations of trap minima relative to 
the APCW with uncertainty below 10 nm). The agreement 
validates the absolute control of our fabrication process (in¬ 
cluding the negligible effect of loss and disorder along the 
APCW), as well as the power of the theoretical tools that we 
have developed lfT4l 1411 l42l . 

We confirm that the variation of ftot in Fig. [^a) is not due 
to the heating of atomic motion during the trap hold time. To 
see this, we adjust N via different MOT loading times and 
measure the decay rate at the shortest hold time (fm = 3 ms), 
as shown by blue points in Fig. [^b) . These observations 
are consistent with those from varying the trap hold time (red 
points in Fig. [^b)), and lead to an almost identical single¬ 
atom decay rate F^qj/Fq = 2.0 ± 0.1 at the shortest loading 
time, corresponding to p/po = 0.16 and N 1. 

The data and our analysis related to Fig. [^strongly support 
the observation of superradiant decay for atoms trapped along 
the APCW. Assuming ftot = fsR -b f^t and fitting ftot lin¬ 
early with N, as shown in Fig. [^b), we find that the supetTa- 
diant rate is given by f sr = ?7 • A^ • Fid with 77 = 0.34 ±0.06. 
The slope 77 is reduced below unity by the random distribution 
of atoms along the x-axis im. 

This observation of superradiant decay is complemented by 
line broadening for steady-state transmission spectra T(A) 
measured at fm = 3 ms with Atm = 5 ms, as show in Fig. 
1^ The measured linewidths f m are significantly broader than 
the free-space width (FWHM) FQ/27r = 4.56 MFIz 1461 . pre¬ 
dominantly due to cooperative atomic coupling to the GM of 
the APCW. We also observe a significant drop in T/Tq at line 
center due to strong atom-photon coupling. Indeed, in Fig. 
(a), we measure T/Tq ~ 0.30 (i.e., a 70% attenuation of the 
GM flux |iTinP) for maximum density po, and T/Tq ~ 0.95 
at the lowest density investigated, p/po « 0.06. 
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FIG. 4: Steady-state transmission spectra r(A) and fitted atomic 
linewidth f m- (a) T{A) with A = 0 corresponding to the free-space 
line center. The three sets of points are measured at relative densities 
p/po ~ 0.12 (black), 0.24 (blue), and 1 (red), where the transmission 
without atoms is Tq. Solid curves are Lorentzian fits to determine 
the linewidth fm- Each point in the spectra is an average over 10 
experiment repetitions, (b) Fitted linewidths (circles) normalized to 
Fo as a function of p/po- The solid line is a linear fit with intercept 

offJ/Vro = 2.1 ±0.1. 

No clear density dependent shift is observed in Fig. I^a), 
in support of our neglect of cooperative energy shifts |iTdd| 
M- The shift in line center for T(A) from A = 0 in free 
space to A = 14 MHz for atoms trapped along the APCW 
is induced by the dipole trap. Furthermore, trapped atoms 
should suffer small inhomogeneous broadening in the spec¬ 
tra shown in Fig. since the FORT shift is small (< 1 MHz) 
for the 6 Pij 2 ,F = 4' excited state, and atoms are well lo¬ 
calized around the trap center due to their low temperature 
T ~ 50pK, corresponding to a small range of light shifts 
< IMHz for atoms in the ground state. 

In Fig. 1^ (b), we plot the linewidths Fm extracted from 
T(A) as a function of p/po- f„i/FQ « 3.4 is largest at 
p/po = 1, and reduces to fm/FQ « 2.1 at p/po = 0.06. 
From linear extrapolation, the single-atom linewidth is esti¬ 
mated to be fm^/FQ = 2.1 ± 0.1. Absent inhomogeneous 
broadening, we expect that fm^ = fiD ± F'. With the cal¬ 
culated F'/Fq « 1.1, the single-atom coupling rate can be 
simply deduced as Fid/Fq « 1.0 ± 0.1. A simple estimate 
of the maximum mean number of atoms then follows from 
Nm = (f„,(pQ) - F')/f ID 2.4 ± 0.4 atoms i?). 

In conclusion, we have used an integrated optical circuit 
with a photonic crystal waveguide to trap and interface atoms 
with guided photons. Superradiance for atoms trapped along 
our APCW has been demonstrated and a peak single-atom 
emission rate into the APCW of Fid/Fq = 1.1 ±0.1 inferred. 
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Our current uniform trap along the APCW is a promising plat¬ 
form to study optomechanical behavior induced by the inter¬ 
play between sizable single-atom reflectivity and large opti¬ 
cal forces (e.g., self organization iiiiiiii). By optimizing the 
power and detuning of an auxiliary guided mode field near the 
air band of the APCW, it should be possible to achieve stable 
atomic trapping and ground state cooling Il48ll49l at trap sites 
centered within the vacuum gap, thereby increasing Fid five¬ 
fold lfT4ll . Opportunities for new physics in the APCW arise 
by fabricating devices with the atomic resonance inside the 
band gap to induce long-range atom-atom interactions ioi- 
1 ^ , thereby enabling investigations of novel quantum trans¬ 
port and many-body phenomena. 
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Supplemental Material: Superradiance for atoms trapped along a photonic crystal waveguide 


DEVICE CHARACTERIZATION 


A schematic of the alligator photonic crystal waveguide (APCW) is illustrated in Fig. ISlja). The waveguide is made from 
200-nm thick stoichiometric SiN with refractive index n = 2.0 IIS2I . The dimensions of the nominal photonic crystals are the 
following: lattice constant a = 371 nm, gap g = 238 nm, width w = 157 nm, and tooth amplitude A = 131 nm, as shown in 
Fig. |STfb). The nominal photonic crystal section consists of A^ceii = 150 unit cells, terminated by 30 tapered cells on each side 
to provide ‘mode-matching’ to and from double nanobeams sections. 

The APCW is characterized by measuring the transmission spectrum Tq{i') without atoms. The resonant structure around 
frequencies I'i displayed in Fig. S2 arises from reflections in the tapered sections at the two ends of the APCW. The free spectral 
range Ai^i = — i'i between resonances decreases as the band edge frequency i^be is approached, which is a signature of an 

increasing group Ug index near i^be, with Ug oc. \/Avi for an ideal structure. 

Our experiment in Ref. im is operated around the frequency Va of the Di: 6 iS'i/ 2 , F = 3 ^ 6 P 1 / 2 , F' = 4 transition in 
atomic Cs, with i'a aligned near i'be by absolute control of the fabrication process at a level of 10“^. Fine tuning for i'a = 1^1 is 
achieved by way of a guided-mode (GM) heating beam with a wavelength of 850 nm and optimum power, typically P > 100/rW. 
In addition, we turn on a strong GM heating beam for 100 ms at the end of each experimental cycle in order to keep the device 
clean by desorbing Cs from the APCW. 

In order to estimate the group index Ug, single-taper reflectivity Rt, and intensity loss we use a model based on the 

transfer matrix formalism for a periodic system to fit the transmission spectrum I1S3I . which we now briefly describe. The 
dispersion relation for the wavevector k{i') near the band edge is approximated by the fitting function | 


k{iy) = ko 


' {eq - + in)) -A2 

-I - A? 


(SI) 


where the wavevector at the band edge is kg = ir/a. Here, fitting parameters are the frequency at the center of the band gap vq, 
the size of the band gap 2Ag, the asymptotic group velocity far from the band edge 2t:v^ jk^, and the loss parameter k. The loss 
parameter k comes from using perturbation theory to add a small imaginary component to the dielectric constant of the material, 
resulting in an imaginary propagation constant that is approximately given by 

Im[fc(i^)] « (S2) 

It provides a convenient way to model losses that scale with inverse group velocity. 

Next we consider the weak cavity formed by the taper reflections Rf. The single-pass phase accumulation (p and single-pass 
power transmission through the cavity are written by 


(p = AceuaRe[fc] and ( = Aceiialm[fc], 


(S3) 



FIG. SI: (a) Schematic of the APCW. An incident field Ein excites the TE-like fundamental mode, and the intensities fro the transmitted 
and reflected field rEin are recorded for device characterization, (b) SEM image of APCW with lattice constant a, gap g, width w, and tooth 
amplitude A. 
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FIG. S2: (a) Measured transmission spectrum Tq(v) for the APCW (black) around the edge of the dielectric band and the model fit (red). The 
dashed lines mark the resonant frequencies Vi from reflections in the taper sections and the solid line marks the band edge frequency t/BE- (b) 
Estimated group index Ug (green) and taper reflection Rt (blue) from the fitted model. For the reference, the transmission spectrum T(,{v) is 
overlaid. At the first resonance marked by the dashed line, the group index is rig « 11, and the taper reflection is Rt ~ 0.48. 


where TVceii is the number of unit cells of the APCW and a is the lattice constant. Then, the transmission through a symmetric 
cavity with mirrors Rt is given by 


^cavity — 


1 

1 + £ + F sin^ [(j)] ’ 


(S4) 


where the coefficient F and loss coefficient L are given by 


e - 2 C (1 - i ?,)2 


and 


F = 


4Ft 

a-Rty 


(S5) 


In order to fit this model to the measured transmission spectrum, first we use the dispersion model (Eq. with K = 0) to fit 
the positions of the cavity resonances to. Second, we fit Eq. ( |S4| i with no loss (C = 0) to the transmitted spectrum by using the 
fitted dispersion model to find (f) and by using a fitting function for F IIS3I . namely 

(i^)-i/2 = Ai(dt./Ag) + A2(dt//Ag)2 + A3(dt./Ag)3, (S6) 


where dv is the distance in frequency from the band edge. Einally, we find the loss parameter C, that makes the on-resonant peak 
heights of the model best match our measurement. 

Eigure [S^ (a) shows the measured transmission spectrum (black curve), overlaid with the model fit (red curve). The fitted 
parameters for the dispersion model are 2Ag = 14.44 THz, v^/vq = 0.60, and vq = 342.8 THz. The fitted parameters for F are 
Ai = 9, A 2 =-48, and A 3 = 128. The fitted loss parameter is k = 1.5 x 10“^. At the first resonance t^i, the model linewidth is 
55 GHz, in reasonable agreement with the measured linewidth of 66 GHz. The fitted dispersion relation is used to estimate the 
group index, and the fitted cavity model is used to estimate Rt and the single pass transmission as shown in Eig. S2 (b). 
At the first resonance, the group index is Ug « 11, the taper reflection is Rt « 0.48, resulting in a peak intensity enhancement 
~ IIS3I . and the single-pass transmission is e~‘^^ » 0.89. Since the propagation loss in the APCW is reasonably 
small, we ignore the loss in our analysis in the following sections. 


FINITE DIFFERENT TIME DOMAIN CALCULATIONS FOR COLLECTIVE COUPLING RATES 


Due to strong coupling to the TE-like GM in the APCW, trapped atoms experience both enhanced atomic decay rates as well 
as collective Lamb shifts. To estimate the size of these effects, we perform EDTD calculations and Eourier analysis as described 
in Ref. IIS4II to obtain the two-point Green’s tensor G(ri, r 2 , to) for the APCW shown in Pig. SI We then evaluate dissipative 
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FIG. S3: (a) Single-atom decay rate (black circles) and excited state level shift (red circles) at = (0, 0, ^i) nm. Vertical dashed 
lines mark the frequencies of the first two guided mode resonances, vi and t/ 2 , near the band gap (frequency range Ua ^ = 335.5 THz 

where Ftot appears constant) that are supported by the finite length of the APCW. Horizontal dashed line indicates F'/Fo = 1.1, estimated 
from the constant F^^j in the band gap region, (b) Dissipative coupling rate Fdd( 2 :) = |F(ra, Va + a:;x)| between two trapped atoms separated 
by X, with their resonant frequencies at either Va = (solid circles) or Va = 336 THz > vbe inside the band gap (open circles), respectively. 
Solid line is an analytical calculation considering actual finite size of the APCW (Fig. [STJ. 


and coherent coupling rates, respectively, as GSHstI 


r(ri,r2) = • Im[G(ri,r2,Wa)] • d 

2 

•/(ri,r2) = -^^^d • Re[Gsc(ri,r2,Wa)] • d, 


(57) 

(58) 


where d is the transition dipole moment, oJa the transition frequency, /tq the vacuum permeability, and h Planck’s constant 
divided by 27r. Here, G^c = G — Gq is the scattering Green’s tensor, in which the vacuum contribution Gq is subtracted from 
the total Green’s tensor G; Im[.] and Re[.] represent imaginary and real parts, respectively. The coupling rate F = r/2 + zJ 
controls collective excitation dynamics of trapped atoms along the APCW. 

We obtain single-atom rates by setting ri = r 2 = at the location of a trapped atom, and evaluate the single-atom total 

, Ua) excited state level shift r, 


decay rate = r(p 


). Figure 


S3 

(T7 


a) shows the calculation 


forr^j = (0,0, zi) nm at the center of the trap shown in Fig. Kb) of ISIII . Here the total decay rate F^^j = Fid + F'(black curve) 
includes the contribution from the GM of interest (Fid), which strongly depends on the atomic resonant frequency i^a = Wo/27r 
and position r^, as well as the coupling rate to all other modes (F'). F' can be estimated from FjQ^(z2a) inside the band gap ( 
^ J^be); F'/Fq « 1.1 remains constant over a broad frequency range. Coupling rate to the TE-like GM, Fid = F^^^ — F', 
can be ob tain ed from this analysis with Fid/Fq = 1.2. 


In Fig. 


S3 


(a), we calculate a small excited state level shift |J^^^|/Fo < 0.4 over a frequency range around vi = 335 THz. 
For our experimental configuration, with « izi, we find | J’^^)(r'o)|/Fo 0. This also suggests that the collective level shift 

for two trapped atoms, |iTdd(a^)| = \ Ta + a:x)| <C (Fq, Fid), is negligible, where x is the atomic separation. Indeed, we 
do not see clear evidence of IV-dependent level shifts in the steady-state transmission spectra shown in Fig. 4(a) of IS II . 

Figure(b) shows Fdd(a:) = |F(rQ,ra + a;x)| for two trapped atoms located at the center of unit cells (x/a G Z) and 
with resonant frequencies at Va = vi or i/a > i'be inside the band gap, where Fj^^ — F' ~ 0. For |a:| > o, Fdd(a:) can be 
used to estimate the dissipative coupling rate between two atoms. When i/a = and |a;|/a > 2, Fdd(a;) slowly drops from 
rdd(O) — F' = 1.2 Fq to smaller values as |x| becomes comparable to the size of the APCW (black circles). This is caused 
by interference with reflections from the tapering regions surrounding the APCW. Solid line in Fig.|S3(b) shows an analytical 
calculation Fdd(a^) = (Fdd(O) — F') cos(7r|a::|/A’effa) that compares to the numerical result, where the fitted effective number 
of cells Nes = 162 ± 9 is larger than Nceu due to the leakage of the fields into the taper regions. Small variations between 
the analytical and numerical calculations are due to residual coupling via other channels. On the other hand, when Va > r'BE 
inside the band gap, Fj^Ca^) quickly drops below O.IFq at |a:|/a > 2. This is expected because, inside the band gap, atoms can 
only cooperatively decay via photonic channels that contribute to F', which are either weakly-coupled or are lost quickly into 
freepace within distances |x| < 2 a. 
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FIG. S4: (a) Normlaized transmission spectmm. Black curve shows the fit to Eq. JS^atfm = 2.5 ms with^ = —0.6 ±0.1, Co = 0.24 ±0.1, 
r = 8.2 ± 0.6 MHz and Aq = 9.3 ± 0.3 MHz. (b) Normalized transmission T/To as a function of the holding time, measured by the on 
resonant guided-mode probe with A = 10.5 MHz. By fitting the measured data to Eq. ( |S9^ with fitted parameters extracted in Eig. |S4| (a), we 
obtain the lifetime of tgm = 28 ± 2 ms (black curve). 


LIFETIME OE TRAPPED ATOMS ALONG THE APCW 


To characterize the lifetime of trapped atoms near the APCW, we measure the normalized transmission T/Tq as a function of 
the measurement time as shown in Fig. 21111 and replotted in Fig. During the lifetime measurement, the frequency 
of the Di transition for the probe field Ein is located between the first and second taper resonances, which leads to the dispersive 
spectrum shown in Fig[^(a). In order to estimate the lifetime of the trap with off resonant cavity, we employ the steady state 
equation IIS8I . 


T/To = (1±0")/ 


1± 


2C(fm) 
1 ± 


+ e- 


2C{U)5„ 

1 + 6^ 


(S9) 


where the normalized detuning from the light shifted resonance Aq is 6^ = 


A-Ar 


, the cooperativity parameter is (7(fm) = 


Co exp(—fm/rcM) with peak cooperatively Co, lifetime tgm, and normalized detuning from taper resonance 9. First, we fit the 
measured spectrum at = 2.5 ms to Eq. ( |S9| l and obtain the fitted parameters, 9 = —0.6 ±0.1, Co = 0.24 ±0.1, F = 8.2 ± 0.6 
MHz and Ao = 9.3 ± 0.3 MHz as shown in black curve in Fig. |S4|(a). Then, to estimate the lifetime tgm, the measured data 
for A = 10.5 MHz in Fig. S4 (b) are fitted to Eq. ( |S9[ ) with fitted parameters from Eig. S4 (a). We obtain the lifetime of 
TGM = 28 ± 2 ms shown in black curve in Eig. |S4|(b). 


MODEL EOR SUPERRADIANCE OF TRAPPED ATOMS 


Our model of superradiance of trapped atoms is obtained by including transfer matrices for atoms in the device model de¬ 
scribed in Section [| OS 31 IS9l4S 111 . Since the first resonance of the taper reflections is aligned to the Di transition for the probe 
Ein, the wavevector of the probe mode is fc = ~ f number of cells Ngs, and the probe field inside the unit cell 

forms a nearly perfect standing wave due to the Bloch-periodic function. In addition, atoms are trapped near the central region 
of the APCW along the x axis (Ax = ±10 /rm). Thus, we ignore the dephasing between atoms and envelope from the taper 
reflections due to the small mismatch of the wavevector Ak = ^ relative to ko = ^ at the band edge. In the following, we 

set the wavevector k = kr, = - and will discuss effects due to the Afc mismatch later. 

The reflection of N atoms randomly distributed at the location Xi with the coupling rate Fid cos^{ kxi) is given by, 

TAf(<5) = where cos^ (kxi). (SIO) 

1 — tCN 1 + 0 ^ 

I 

where the single-atom fractional coupling rate is = Tid/F' and normalized detuning is (5 = 2A/r' with F'/Fg « 1.1 
from the numerical simulation in Section lEl. The temporal profile of superradiance from N atoms is obtained by Fourier 
transforming r]s[{5) to yield rAr(f), and taking the convolution of rpf{t) with a gaussian pulse of the half width tr ^ 5 ns for the 
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FIG. S5: (a) Temporal profiles of atom ic em ission into the GM at tm = 63 ms with and without atoms, shown in cyan and black circles, 
respectively. The blue curve shows Eq. IS12 1 fitted for a single atom with Tid/T' = l.OiO.l. (b) Temporal profiles of atomic emission into 
the GM at fni = 3 ms with and without atoms, shown in pink and black circles, respectively. The red curve shows Eq. ( |S13| > fitted to yield 
with N = 2.6 ± 0.3. The background level of (a) is higher than (b) due to the drift of the intensity modulator during the 5 times longer data 
accumulation time. 


excitation pulse Furthermore, the temporal profile at f > 2cr can be approximated by 


'{t) = 


' N 




_^ V 2 

riD E cos^ (kxi) j ■ exp 


- ( E + r' j f 


(Sll) 


Considering the assumed random locations of atoms with uniform probability density in the unit cells along x, the spatially 
averaged temporal profile is obtained by integrating Eq. (Sll[) along x, yielding 


lN{t) = 7^e • Jo (7f) 


N-2 


- ^ - -Io {-it) - \^— + — j Io {it) h (jt) + - - - -Ii (jt) 


iS12) 


where Ik(z) is a modified Bessel function of the first kind and 7 = riD/2. In addition, the number of trapped atoms along the 
APCW is drawn from a Poisson distribution p(JV, N) with mean number of atoms N. The total decay curve then becomes 


1 -tot{t) = Co'^^p{N,N) •l7v(f) +Ibg, 


(S13) 


N 


where cq is a constant. Here, the background intensity Ibg is measured separately without atoms and is given by the black 
circles in Fig. S5 As shown in Fig. 3 flUl, the total decay rate asymptotes to r^Qj/Fo = 2.0 ± 0.1 at the longer hold times, 
which suggests that atomic decay at fm = 63 ms mostly originates from a single atom for mean atom number N 1. Thus, we 


fit the decay curve at f„i = 63 ms toEq. (S12lwith7V = 1, and obtain ^0 = Tio/r^ = 1.0 ± 0.1 shown in Fig. S5 (a). Then, by 
using the fitted ^ 0 ^ the shortest hold time data {tm = 3 ms) is reasonably well fitted to Eq. (S131 with N — 2.6 ± 0.3, as shown 
inFig.[sg(b). 

We also numerically estimate the contribution of the envelope from the taper reflections and dephasing between atoms 
along the APCW. We employ the transfer matrix model with k = ko — Ak, which includes the coupling rate riD(a;) ~ 
cos^(fcoa;) cos^{Akx) along the x-axis of the APCW and propagation phase Ak ■ 6x between atoms separated by Sx. Here, 
denotes the peak coupling rate of both the unit cell and the envelope from taper reflections. The position of the atoms is 
generated from a normal distribution with = 10 pm at the temperate of 50/iK. Then, we numerically generate the decay 
curve and extract the total decay rate, suggesting that F id extracted from Eq. ( |S12| i with = 1 is underestimated by ^ 10% 
and IV from the fits to Eq. (S13lby ^ 15%. Note that we do not incorporate these corrections in our estimation of Fid and IV 


in Ref. Ell, since the temperature of the atoms trapped along the APCW could be different from the measured temperature in 
free space. Indeed, the calculated trap potential combined with Casimir-Polder potential suggests that atoms trapped along the 
APCW could be much colder (< 20 pK) than the measured temperature for the free-space FORT (~ 50 pK) due to the smaller 
trap depth near the APCW, leading to a smaller correction of Fid and N due to tighter localization of the atoms around the 
center of the APCW. The distribution of atoms along the APCW is being investigated in more detail. 

To support ou r assu mption of TV-dependent superradiance, f tot = f sr + f id with f sr = p • TV • Fid, we generate the decay 
curve from Eq. (S13 1 with Fid/Fq = 1.0 and various TV, and extract ftot by fitting to an exponential. The dashed curve in Fig. 
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N 

FIG. S6: Fitted total decay rate Ftot normalized by Fq as a function of the mean number of trapped atoms N. The dashed curve shows the 
calculated Ftot from the model, overlaid with results measured for various hold times (red circles) and for loading times (blue circles). A 
linear fit to the combined data (solid black line) gives Fsr = ry • iV • Fid with r} = 0.34 ± 0.06. The inset shows a log-log plot of the curve 
generated from the model. 


S6 shows the calculated ftot, overlaid with measured hold time (red circles) and loading time (blue circles) dependence and the 
linear fit (solid black line). Although the dashed curve generated from the model deviates from the linear dependence at iV < 1, 
the linear fit captures the N dependence reasonably well. As clearly seen in the inset of Fig. S6 the nonlinear dependence on 
TV at (V <C 1 is due to the “conditional” character of decay rate measurements, meaning that the decay curve consists mostly of 
fluorescence from a single atom, despite N 1. Due to the negligible background counts in our measurements, single detection 
events at (V ^ 1 herald the presence of single atoms. 

The linear fit to the combined data sets gives f sr = ■ Fid with rj = 0.34 ± 0.06, consistent with the model for N > 0.7. 

A qualitative understanding of this value of rj is the following. Due to the random distribution o f ato ms along the APCW, the 
intensity of atomic emission into the GM is spatially modulated by cos^(A:a:) as shown in Eq. (Sill, meaning that both GM 
excitation of atoms and emission into the GM are proportional to cos^{kx), resulting in cos^(fcx) dependence. With the spatial 
averaging along the x-axis of the APCW, the superradiant decay rate is then reduced by a factor of roughly p ^ 3/8 (i.e., the 
average of cos^(A:x) over a unit cell). 
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